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We implement the discrete-time quantum walk model using the continuous-time evolution of the Hamiltonian 
that includes both the shift and the coin generators. Based on the Trotter-Suzuki first-order approximation, we 
consider an optimization problem in which the Hellinger distance between the walker probability distributions 
resulted from the evolution of such Hamiltonian and the quantum walk dynamics is minimized. The phase space 
implementation of the quantum walk is considered where the walker state is encoded on the coherent state of a 
resonator and the coin on the two-level state of a qubit. In this approach, no mechanism for switching between 
the coin and the shift operators is included. We show the Hellinger distance is bounded for large number of time 
steps. The distance is small when we deviate from the standard quantum walk model, namely, when the walker 
is allowed to move in between the sites. In simulating the standard quantum walk model, the distance is large 
but bounded by 26% for a relevant number of time steps. Even so, the system evolution shows the essential 
characteristics of the standard quantum walk dynamics, namely, the ballistic evolution of the probability distri¬ 
bution and the linear growth of the corresponding standard deviation. Moreover, the entanglement generated in 
this approach is approximately the same as the entanglement generated in the standard quantum walk dynamics. 
Finally, the system dynamics under the influence of the decoherence also shows the similar quantum-to-classical 
transition as in the standard quantum walk dynamics. 

PACS numbers: 03.67.Ac,05.40.Fb,42.50.-p,02.30.Yy 


I. INTRODUCTION 

Discrete-time quantum walks, the counterpart of the clas¬ 
sical random walks, have recently attracted considerable at¬ 
tention in designing quantum search algorithms that outper¬ 
form the similar classical algorithms H]. It has been also 
shown that the discrete-time quantum walks are universal for 
quantum computation El. Many experimental setups, there¬ 
fore, have been proposed to implement quantum walks and 
several implementations have also been reported 0. Ex¬ 
amples of systems explored include micromasers 0], ion 
traps (Ml, cavity quantum electrodynamics IHl, ensembles 
of nitrogen-vacancy centers in diamond 121 and optomechan¬ 
ical systems Go). 

In order to implement the discrete-time quantum walk 
model, two separate unitary operators are required, namely, 
the coin and the shift operators. The corresponding Hamilto¬ 
nians do not commute hence the quantum walk dynamics is 
generated by alternating their actions. This is, however, tech¬ 
nically difficult to implement depending on the physical sys¬ 
tem at hand. The discrete-time quantum walks may also be re¬ 
alized by using an always-on shift Hamiltonian but applying a 
piecewise coin Hamiltonian from time to time. A sequence of 
pulses is then required to implement the coin operator. The ap¬ 
plication of such pulses, however, generally introduces noise 
and disturbs the walker’s position. 

In this way, the experimental realizations with the 
continuous-time quantum walk model im, that involves no 
coin operator, are usually more convenient and preferred. 
However, there is no known search algorithm based on the 
continuous-time model that outperforms any classical algo¬ 
rithm, when the dimension of the data structure lattice is less 
than four Ha. For the two- or three-dimensional lattices, effi¬ 
cient algorithms exist just whenever the lattice spectrum pos¬ 
sesses Dirac points EMnl. A strategy to overcome the prob¬ 


lems of the implementation of the discrete-time quantum walk 
model, would be, therefore, to keep the discrete-time quantum 
walk model but implementing the walk with the continuous¬ 
time evolution of the system. The same strategy may be ap¬ 
plied to coinless quantum walks, for example to the one de¬ 
scribed in ifT^ . 

In this paper, we employ the Trotter-Suzuki first-order ap¬ 
proximation lfT9l I 20 II . exploring the situations in which the 
discrete-time quantum walk model can be efficiently simu¬ 
lated by the continuous-time evolution of the total Hamilto¬ 
nian that includes both the coin and the shift generators. In 
fact, the continuous-time evolution of the total Hamiltonian 
when considered at some specified time steps approximately 
produces the discrete-time quantum dynamics, with a variable 
amount of error. Having realized the coin and the shift oper¬ 
ators in a single time step, we call this approach the discrete¬ 
time quantum walk with simultaneous coin and shift opera¬ 
tors. In such implementation, the shift and the coin genera¬ 
tors in the Hamiltonian are always on. Even so, the effects 
of the non commutativity of those generators does not destroy 
the important features of the discrete-time quantum walk dy¬ 
namics. Those features include the ballistic evolution of the 
walker probability distribution, the linear growth of the cor¬ 
responding standard deviation, the amount of the generated 
coin-walker entanglement and the quantum-to-classical tran¬ 
sition of the probability distribution in the presence of the de¬ 
coherence. The significance of such implementation in de¬ 
signing efficient quantum algorithms, however, needs to be 
addressed elsewhere. 

We specifically consider the quantum walk on a circle in 
the phase space with a two-sided coin, given its simplicity 
and potentiality for implementation on a variety of physical 
systems EHini, but our findings extends to other geometries 
as well. After a brief review of the standard implementation of 
the discrete-time quantum walk on a circle in Sec. we show 
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the total Hamiltonian including both the shift and the coin 
parts can approximately simulate the quantum walk dynamics 
up to a certain fidelity in Sec. Ill By modifying the Trotter- 
Suzuki approximation we optimize the system frequencies to 
minimize the approximation error. Based on the walker jump 
at each time step we introduce a parameter that specifies the 
level of the error. We finally consider more closely other char¬ 
acteristics of both the continuous-time evolution and the stan¬ 
dard discrete-time quantum walk dynamics, such as the stan¬ 
dard deviation, the coin-walker entanglement, as well as the 
effect of a decoherent coin in Sec. IV Finally in Sec. [V|we 
present our conclusions. 


II. THE STANDARD QUANTUM WALK ON A CIRCLE 

The dynamics of the discrete-time quantum walk is de¬ 
scribed by 


\i:{t + l))=S{I^C)\m), ( 1 ) 

where S and C are the shift and the coin operators, respec¬ 
tively, and I is the walker space identity. At each time step, 
the state of the quantum walk \ip{t)) is a superposition of the 
walker-coin compound states. For a walk on a circle with d 
sites the generic state of the quantum walk at time t is given 
by 


d-1 

\^m}, ( 2 ) 

m—0 

where ^o,i the probability amplitudes in terms of the site 
basis {\ipm = 27rm/fi), TO = 0 ... d — l} and the coin basis 
{|0), |1)}. The interesting quantity here is the walker proba¬ 
bility distribution over the sites that is obtained by tracing out 
the coin space. Unlike the behavior of the classical random 
walk, the dynamics of the probability distribution here is bal¬ 
listic and the corresponding standard deviation grows linearly 
with the time steps m. 

In the phase space realization of the quantum walk, the 
walker is encoded on the coherent state of a resonator field 
and the coin on the two-level state of a qubit which is coupled 
to the resonator. The required Hamiltonian is 

7i = -huiax — hga) aaz, (3) 

where w is the qubit transition frequency, g is the qubit- 
resonator coupling constant, (a) is the resonator creation 
(annihilation) operators and and are the Pauli matrices. 

Suppose the initial state of the system is given by 
|'0o) = |a)k) where |a) is the coherent state of the res¬ 
onator and \q) is the state of the qubit. Suppose, furthermore, 
there is a mechanism that allows to switch between the first 
(the coin) and the second (the shift) terms in Hamiltonian ([^, 
alternately, during the evolution time t. The challenge in 
the experimental realizations of the discrete-time quantum 
walk model is, actually, to effectively establish such an “al¬ 
ternating” Hamiltonian. Nevertheless, consider a uniform 


time stepping of duration tjn = t with which the Hamil¬ 
tonian is switched. By choosing the resonator frequency such 
that ojT = 7r/2, the first term in Hamiltonian © generates 
a Hadamard-like transformation (7r/2-rotation around the x 
axis) for the qubit, thus the coin operator C is implemented. 
The second term generates a rotation for the coherent state in 
the phase space by the angle igr whose sign depends on the 
state of the qubit. In this way, the shift operator S is realized. 
Therefore, after each period 2r a single step of the quantum 
walk has been completed. 

To explore the system dynamics, a finite dimensional phase 
space of size d is required which is obtained by truncating the 
infinite dimensional Fock space. The number operator a^o 
is expected to be conjugate to the phase operator. Therefore, 
the basis for the finite-dimensional phase space is constructed 
by calculating the discrete Fourier transform of the number 
states in the truncated Fock space. The initial coefficients 
in Eq. are determined by calculating the dis¬ 
crete Fourier transform of the initial truncated coherent state 
and knowing the initial coin state. Note that the initial coher¬ 
ent state is a non-local state in the phase space with a distribu¬ 
tion close to the Gaussian. The dynamics of the system with 
the “alternating” Hamiltonian is then obtained, using Eq. Q- 
The shift operator in terms of the phase space basis is given by 
where is the Eourier transform operator. Choosing 
the coupling constant such that gr — 2Tr/d forces the walker 
to hop between the sites on the circle. The phase probability 
distribution of the walker, at any time step, is obtained by con¬ 
sidering the diagonal elements of the density operator of the 
system after tracing out the coin space. Experimentally, op¬ 
tical homodyne tomography can be applied to reconstruct the 
state of the walker, if it is associated to an optical field, and 
ultimately used to reconstruct the phase probability distribu¬ 
tion ETll . 

The truncated dimension d which is also the total num¬ 
ber of sites for walking on the circle is highly restricted in 
Ref i). There, it is demanded the sites are well separated 
in order to have approximately complete overlap between the 
initial coherent state and a single phase state, hence to reach 
a local initial walker state. Therefor, by fitting 27r|a| distin¬ 
guishable coherent state on a circle of radius |a|, the strict 
range |ap + |a| < d < 27r|a| for the truncated dimension 
is established. The lower bound is required for the coherent 
state has reasonable support on the truncated Eock (or phase) 
space i). Such restrictions leads to the maximum number 
of sites dmax = 33 for |a|max = 27r — 1. However, treat¬ 
ing the coherent state as a non-local initial state in the phase 
space makes the upper bound for the dimension d unneces¬ 
sary, as far as the quantum walk dynamics is concerned. But, 
if just the distinguishability of the coherent states on the cir¬ 
cle is required a wider interval can be considered. Actually, 
approximately 47r|Q;| distinguishable coherent state can be fit¬ 
ted on a circle with radius |a| Go). That can be justified 
by considering the geometrically obtained uncertainty relation 
A(pAn > 0.506 for |a| > 1. Here, Aip is the phase uncer¬ 
tainty and An = |a| is the field excitation number uncertainty. 
In this paper, we consider the wider interval which leads to the 
maximum number of sites dmax = 145 for |a|niax = 47r — 1. 
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It should be mentioned that for a given |a| choosing the largest 
allowed value of d does not always lead to the best quantum 
walk dynamics. Such cases, however, can be cured by chang¬ 
ing the dimension. 


III. SIMULATION WITH SIMULTANEOUS COIN AND 
SHIFT OPERATORS 

The above scenario is considered as the ideal way for im¬ 
plementing the discrete-time quantum walk model. The basic 
problem in implementing such method is the technical diffi¬ 
culties in realizing a mechanism for switching between the 
shift and the coin terms in Hamiltonian ([^. In this section, 
we analyze the simulation of the discrete-time quantum walk 
dynamics using simultaneous coin and shift operators. That 
can be realized by continuous evolution of the system without 
needing any switching mechanism. 

The time evolution of Hamiltonian ( [^ ca n be decomposed 
using the Trotter-Suzuki approximation i 1 9l l20ll 

^ _ g-i'Ht/S ^iga^aa^t/n 

where the approximation error scales as guit'^/n that can be 
made arbitrarily small by increasing n. The approximated 
evolution in Eq. 0 corresponds to a n-step discrete-time 
quantum walk dynamics. At each time step t/n = t, the state 
of the qubit is transformed by the rotation TZx{ujt) about the 
X axis, hence the coin operator is realized. Depending on the 
state of the qubit the coherent state rotates by the angle ±pr 
accordingly, thus the shift operator is also realized (see IH). 
In this way, the operator U in Eq. 0 which is correspond¬ 
ing to the continuous-time evolution of the system seems to 
be useful for simulating the discrete-time quantum walk dy¬ 
namics. Actually, the continuous-time evolution considered 
at time steps r mimics a discrete-time quantum walk dynam¬ 
ics. Note that while in the case of “alternating” Hamiltonian 
a single step of the quantum walk is completed in period 2t, 
in the case of Trotterized Hamiltonian it is completed in r. 
Therefore, associating the latter case with the quantum walk 
with simultaneous coin and shift operators seems reasonable. 

By setting the system frequencies and the time step such 
that gr = 2^/d and ujt — 7r/2, the approximated evolu¬ 
tion in Eq. 0 generates a quantum walk dynamics with the 
Hadamard-like coin. Such dynamics, however, does not cor¬ 
respond to the exact evolution of the system in Eq. Q. In fact 
with the mentioned setup, the walker phase probability distri¬ 
butions resulted from the exact and the approximated evolu¬ 
tions will be very different after a few time steps. We mea¬ 
sure the difference between the probability distributions by 
the Hellinger distance ll22]| 

i?(/^ex,Papp)= (5) 

where Pex and Tkpp are the walker phase probability dis¬ 
tributions at a given step of the evolution, corresponding to 
the exact and the approximate dynamics, respectively. The 


norm || .||2 is the Euclidean norm. We associate the distance in 
Eq. (|^ with the approximation error in Eq. Q. 

The approximation error in Eq. 0^ can be decreased by 
decreasing the time step r. But, decreasing t will specifically 
decrease the angular rotation of the qubit ujt, hence a coin 
operator with ujt < 7r/2 is realized. Specially, when a very 
small upper bound error is allowed, r is required to be very 
small and therefore the coin gets very close to the identity. 
With such a coin, the dynamics is almost deterministic and 
the initial probability distribution is just displaced at each time 
step without deformation. 

It is possible to decrease the error and still generate the op¬ 
timal desired Hadamard coin by decreasing g but keeping ui 
and T fixed. The error, however, will not go to zero just by de¬ 
creasing g. Very small p’s, make the dynamics very slow with¬ 
out really improving the approximation when the total evolu¬ 
tion time is of order r. This is the time order within which 
the ballistic spread of the probability distribution is occurred. 
In fact, the condition ujt = 7r/2 for obtaining the Hadamard 
coin causes the Trotter-Suzuki approximation no longer pro¬ 
duces accurate results in time scales of r. The reason can 
be explained by considering the Zassenhaus expansion of the 
system dynamics f23\ . In the Zassenhaus expansion of the ex¬ 
act evolution in Eq. 0- which is actually an infinite product 
of exponential operators, exponents like (a;T)^(pr)* appears 
where k and I are positive integer numbers. Therefore, with a 
Hadamard coin, it is not possible to obtain a good approxima¬ 
tion by keeping just the first two exponentials in the expansion 
no matter how small g is. 

However, in the above case where g is decreased and the 
other parameters are kept fixed, the difference between the 
phase probability distributions resulted from the exact and the 
approximated dynamics in Eq. Q seems to be mainly due 
to the displacement of the distributions. To improve the ap¬ 
proximation we therefore slightly modify the Trotter-Suzuki 
approximation in a single time step 

^i(gTa''aa,, - u>Taj,/2) ^ ^(2 t! jrd)o'aa^ g-i(-7r/4)tTj, 

where the dimensionless frequencies pr and ujt need to be 
optimized such that the approximation error is minimum. The 
parameter r is a positive number that is used for controlling 
the angular step of the simulated quantum walk. The idea here 
is to adjust the system frequencies such that the continuous¬ 
time evolution of the system coincide the desired discrete-time 
quantum walk dynamics at time steps r, as mush as possible. 
The dimensionless frequencies are optimized by minimizing 
the average Hellinger distance in Eq. 0 calculated for large 
number of time steps. To be close enough to the global op¬ 
timal solution, the maximization process is carried over 100 
random initial guesses for the frequencies and then the solu¬ 
tion with the smallest distance is selected. 

Figure shows the Hellinger distance between the phase 
probability distributions resulted from the exact and the ap¬ 
proximated evolutions in Eq. (j^, for different values of r. 
Two cases with the phase space dimensions d = 31,125 and 
corresponding average field excitation numbers |a| = 5,10 
have been considered here. In both cases the initial state of 
the qubit was set to the ground state jg) = |0). Moreover, 
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FIG. 1. (Color Online) The Hellinger distance between the two prob¬ 
ability distributions resulted from the exact and the approximated 
evolutions in Eq. in terms of the time steps. The distances are 
given for five values of r in the case d = 31 and a = —5 in panel 
(a) and the case d = 125 and a = —10 in panel (b). In both cases 
the initial state of the qubit is set to the ground state. 


the total number of steps at each r is set such that the initial 
distribution is split into two peaks with the angular distance 
of TT radian. The diagrams show that by increasing r the ap¬ 
proximation error over the time steps decreases. Increasing 
r can be associated with decreasing g while keeping fixed w 
and r, that improves the approximation as described before. 

As previously mentioned, the phase space dimension d cor¬ 
responds to the total number of sites on the circle and spec¬ 
ifies the walker Hilbert space dimension. Comparing panels 
(a) and (b) in Fig. shows that the effect of increasing the 
dimension d is just to diminish the oscillations in the error 
plots. 

According to Fig.[T] for the case r = 10 (green/lowest plots 
in both panels) the approximation error is almost stable and 
remains below 4% for a large number of time steps. However, 
since the angular distance between the sites is 27r/(i, choos¬ 
ing r = 10 causes the walker is conditionally shifted by 1/10 
of the distance between the sites. So, in spite of the good 
match between the exact and the approximated evolutions in 
this case, the approximated evolution does not represent the 
standard quantum walk evolution. In the standard model the 
walker is restricted to jump between the neighbor sites at each 
time step. Nevertheless, the standard deviation of the phase 
probability distribution in both exact and approximated dy¬ 
namics is still linearly proportional to the number of steps but 
the corresponding slopes drop at least one order lower, com¬ 
paring the standard quantum walk case. 

For the case r = 1 (magenta/highest plots in both panels) in 
Fig.[T] the approximated dynamics in Eq. corresponds to 
the standard quantum walk evolution. In this case, however, 
the error is larger but remains below 26% until the initial dis¬ 
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FIG. 2. (Color Online) The phase probability distributions for the 
exact (solid line) and the approximated (dashed line) evolutions in 
Eq. ([§ with r = 1, when the peaks are separated by tt radian. Panels 
(a) and (b) correspond to the cases d = 31 (a = —5) and d = 125 
(a = —10), respectively. In both cases the initial state of the qubit 
is set to the ground state. 


tributions are split into two peaks with the angular distance of 
TT radian. To get an idea of how the two dynamics are different 
at such error the corresponding phase probability distributions 
are presented in the next figure. 

Figure shows the phase probability distributions for the 
exact (solid line) and the approximated (dashed line) evolu¬ 
tions corresponding to the case r = 1 in Eq. (|^, for two 
different dimensions d = 31,125. The plots show the prob¬ 
ability distributions for the time step when the two peaks are 
separated by the angular distance of tt radian. It can be seen 
that, the probability distributions are qualitatively similar, rep¬ 
resenting two major peaks. The exact dynamics, therefore, 
represents one of the important features of the quantum walks, 
the ballistic spread of the walker probability distribution. A 
special difference between the two plots corresponds to the 
intermediate interference between the two peaks related to the 
approximated evolution (the standard quantum walk dynam¬ 
ics). Eor the exact evolution the probability distribution is al¬ 
most flat between the two peaks. 

Actually, it is not possible to reproduce the intermediate in¬ 
terference with just playing with the system frequencies. It is, 
nevertheless, tempting to think of including more terms in the 
Hamiltonian of the exact evolution in Eq. to improve the 
approximation. Specially, higher order terms from the Baker- 
Campbell-Hausdorff (BCH) formula ll23l may seem suitable 
in this respect. However, according to the convergence theo¬ 
rem proved in Ref. Il23l the BCH formula does not converge 
because of the large argument lot = 7r/2 of the Hadamard 
coin. Therefore, the higher order terms from BCH formula 
does not lead to an improved approximation. Even so, we 
have numerically verified that adding a term proportional to 
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a^aax can improve the approximation to some degree. It may 
be stressed that this term is different from the second order 
term that comes from the BCH formula and is proportional to 
a^aay. We are not going to consider such term in this paper. 

We now investigate in more details the case r = 1 in 
Eq. (§, which corresponds to the standard quantum walk dy¬ 
namics. Panels (a) and (b) in Fig. show the Hellinger dis¬ 
tance between the exact and the approximated probability dis¬ 
tributions for large number of time steps, in the two cases 
d = 31,125. The plots in this panels are actually the contin¬ 
uations of the plots in Fig.[2 for the case r = 1. The system 
frequencies for producing all the plots have been optimized 
by minimizing the average Hellinger distance over the first 12 
time steps for the dimension 31, and the first 46 steps for the 
dimension 125. However, including more time steps in the 
optimization process will lead to a more balanced error plots, 
specially in the case d = 31 presented in panel (a). 

For time steps less than the total number of sites on the cir¬ 
cle, the exact dynamics shows no interference (see Fig. |^. 
However, the bounded error in Fig. indicates that as time 
step goes beyond the total number of sites on the circle the in¬ 
terference is generated. We have verified (not showed here) 
the two major peaks in Fig. (solid lines) really interfere 
when pass through each other. 

The standard deviations of the probability distributions re¬ 
lated to Eq. (j^ are showed in panels (c) and (d) in Fig. for 
the two cases c? = 31,125. The solid/dashed line corresponds 
to the exact/approximated evolution. After the initial distribu¬ 
tion is split into two peaks with the angular distance around 
27r radian, the peaks will move toward each other hence the 
standard deviation starts to decrease. The oscillations of the 
standard deviation in the hgure indicates the quantum walk 
spreads around the circle for several times. There is a good 
match between the two plots at least for time steps less than 
the available sites on the circle. The exact evolution, therefore 
represents another important feature of the quantum walks, 
the linear growth of the standard deviation of the walker prob¬ 
ability distribution. 

Finally, the entanglement between the qubit and the res¬ 
onator is given by the negativity in panels (g) and (h) in Fig.[^ 
The thick/thin line corresponds to the exact/approximated 
evolution in Eq. (j^. The figure shows there is no important 
difference between the entanglement that is generated during 
the evolution in both the exact and the approximated evolu¬ 
tions. The entanglement dynamics is just smoother for the 
exact evolution. 

Figures and 1^ therefore, suggest the discrete-time quan¬ 
tum walk with simultaneous coin and shift operators is 
promising in simulating the standard discrete-time quantum 
walk model. Despite some differences, the implementation 
with simultaneous coin and shift operators have the important 
features of the standard quantum walk, namely, the ballistic 
spread of the walker probability distribution and the linear de¬ 
pendency of the corresponding standard deviation on the time 
steps. Moreover, the amount of the generated entanglement 
between the coin and the walker is almost the same level as 
the one generated in the standard quantum walk model. 

Inspecting the optimized system frequencies in Eq. •S' re¬ 


veals that the ratio 

^ UJT/gr ^ ^ 

(7r/2) / (27r/rd) rdg 

approximately amounts to 10, in the case of r = 1. As the 
optimized value of gr is not very different from 2Tr/d the 
optimized value of lot is about Cn 12. This is true for both 
dimensions d = 31 ,125 as well as the other relevant dimen¬ 
sions. Therefore, the system frequencies can be adjusted very 
easily to obtain an approximate quantum walk dynamics with 
the exact evolution in Eq. 'll])- We have verified numerically 
the ratio C depends on the coin that is used in the quantum 
walk dynamics, namely the qubit rotation TZx{u!t) around the 
X axis. The value of C changes in the direction of changing 
the angle of rotation. Decreasing the angle from tt/2 (chang¬ 
ing the coin toward the identity) decreases C, and in the con¬ 
trary, increasing the angle from ■k 12 (changing the coin toward 
the completely biased coin) increases C. Actually, the differ¬ 
ent scales of the system frequencies allows to develop the de¬ 
sired qubit rotation. In this way, it is possible to implement 
any desired coin. 

IV. DECOHERENT COIN 

As mentioned before, the probability distribution of the 
quantum walk spreads ballistically and the corresponding 
standard deviation evolves linearly in terms of the time steps. 
In the presence of the decoherence, however, the quantum 
walk makes a transition to a random walk in which the spread 
of the probability distribution is Gaussian and the standard de¬ 
viation evolves with the square root of the time steps. 

Simulating the effect of the dephasing channel Il24l on the 
two-sided coin (the qubit) for both the exact and the approx¬ 
imated dynamics in Eq. (|^ shows qualitatively the same be¬ 
haviors for the case r = 1. For a given dephasing channel, 
the Hellinger distance between the corresponding probabil¬ 
ity distributions decreases when r increases. Moreover, the 
Hellinger distance is decreasing in terms of the time steps for 
nonzero decoherence. Actually, in the classical limit both the 
exact and the approximated dynamics become identical. 

V. CONCLUSIONS 

In this paper, we have analyzed the simulation of the 
discrete-time quantum walk dynamics, using the continuous¬ 
time evolution of a coupled resonator-qubit system. The co¬ 
herent state of the resonator is served as the walker and the 
state of the qubit as the two-sided coin. Instead of successive 
applications of the coin and the shift operators, we apply a sin¬ 
gle operator which is generated by the evolution of the Hamil¬ 
tonian that includes both the coin and the shift generators. Ac¬ 
tually, the continuous-time evolution of the system considered 
at the specihed time steps imitate the discrete-time quantum 
walk dynamics. We have called this method the discrete-time 
quantum walk with simultaneous coin and shift operators. In 
this approach, no mechanism for switching between the coin 
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FIG. 3. (Color Online) Different characteristics of the exact and the approximated evolutions in Eq. for the case r = 1. The left and 
the right panels correspond to the cases d — 31 {a = —5) and d = 125 (a = —10), respectively. Panels (a) and (b) show the Hellinger 
distance between the corresponding probability distributions for large number of time steps. Panels (c) and (d) show the standard deviations 
of the probability distributions over the time steps for the exact (solid line) and the approximated (dashed line) evolutions. The negativity 
entanglement between the qubit and the resonator for the exact/approximated evolutions is showed by the thick/thin line in panels (e) and (f). 


and the shift operators is included. In fact, the corresponding 
operators are implemented in a common single time step. 

Based on the Trotter-Suzuki approximation, we wrote an 
optimization problem whose global solutions (the system 
frequencies) made the approximation error minimum. The 
approximation error is calculated by the Hellinger distance 
between the walker probability distributions resulted from 
the system dynamics and the quantum walk dynamics. We 
showed the error is bounded for large number of time steps. 
The error is small when we deviate from the standard quan¬ 
tum walk model, namely, when the walker is allowed to move 
in between the sites. For the standard quantum walk model, 
the error is large and bounded by 26% for a relevant number 
of time steps. 

In spite of that error, the system evolution with simultane¬ 
ous coin and shift operators shows the essential characteristics 
of the standard quantum walk dynamics. The probability dis¬ 
tribution spreads ballistically and the corresponding standard 
deviation grows linearly in time. The entanglement generated 
in this method is also approximately the same as the entan¬ 


glement generated in the standard quantum walk dynamics. 
Moreover, the system dynamics under the influence of the de¬ 
coherence shows the similar quantum-to-classical transition 
as in the standard quantum walk dynamics. 


It is expected the discrete-time quantum walk model with 
simultaneous coin and shift operators to be useful in imple¬ 
menting quantum algorithms that outperforms classical ones. 
It has been proved the discrete-time quantum walks are uni¬ 
versal for quantum computing Q. Actually, the quantum cir¬ 
cuit model can be transformed into a graph on which the quan¬ 
tum walk propagates and implements the quantum algorithm. 
We expect the similar property for the quantum walk with si¬ 
multaneous coin and shift operators. However, at least a two- 
dimensional quantum walk is required to implement universal 
quantum computation. The extension of the above model to 
more than one dimension is feasible, in principle. It is there¬ 
fore interesting to check in what sense would the above model 
be universal for quantum computing. 
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